Images of a Bose-Einstein condensate: diagonal dynamical Bogoliubov vacuum 
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Evolution of a Bose-Einstein condensate subject to a 
time-dependent external perturbation can be described by 
a time-dependent Bogoliubov theory: a condensate initially 
in its ground state Bogoliubov vacuum evolves into a time- 
dependent excited state which can be formally written as 
a time-dependent Bogoliubov vacuum annihilated by time- 
dependent quasiparticle annihilation operators. We prove 
that any Bogoliubov vacuum can be brought to a diagonal 
form in a time-dependent orthonormal basis. This diagonal 
form is tailored for simulations of quantum measurements on 
an excited condensate. As an example we work out phase im- 
printing of a dark soliton followed by a density measurement. 

Quantum measurements on Bose-condensed systems 
can give quite unexpected results. For example, in the 
classic paper by Javanainen and Yoo [1] a density mea- 
surement on a Fock state \N/2, N/2) with TV particles 
equally divided between two counter-propagating plane 
waves e^" reveals an interference pattern p{x\(p) ~ 
cos^(a; — (fi) with a phase ip chosen randomly in every 
realization of the experiment. The Fock state has a uni- 
form single particle density distribution, but its measure- 
ment unexpectedly reveals interference between the two 
counter-propagating condensates. The Fock state is a 
quantum superposition over A^-particle condensates with 
different relative phases ip in their wave functions [2], 
\N/2,N/2) Jdifi \N : e+^t^-"^) -he-^^^"^)), but every 
single realization of the experiment reveals such a den- 
sity distribution as if the state before the density mea- 
surement were one of the condensates : e+*^^~'^-' -I- 
^-i{x-ip)^ with a randomly chosen phase ip. This ef- 
fect is best explained [2] when the density measurement, 
which is a destructive measurement of all particle po- 
sitions at the same time, is replaced by an equivalent 
sequential measurement of one position after another. 
With an increasing number n of measured positions a 
quantum state of the remaining N — n particles grad- 
ually "collapses" from the initial uniform superposition 
over all phases to a state with a more and more local- 
ized phase (p. For a large N a measurement of only a 
small fraction <C 1 of all particles practically collapses 
the state of remaining N — n particles to a condensate 
with phase p. The lesson from this instructive exam- 
ple [1,2] is that when we want to predict possible out- 
comes of measurements on a Bose-condensed state, then 
it is helpful to rewrite the state as a superposition over 
condensates J Vcj) \N : 0) with different conden- 

sate wave functions (j){r) in a hope that P{(f>) ^ 
will give an approximate probability for different density 
measurement outcomes p{r\<j>) — N\(j>{r)\'^. However, as 



the set of all condensates \N : ip) is non-orthogonal and 
overcomplete, the "amplitude" "0(0) is not unique but 
depends on coordinates used to parameterize the space 
of condensate wave functions In general a special 

care should be taken to choose the right parameteriza- 
tion where the simple "Born rule" -P(0) ^ |V'(0)p gives 
the correct probability. 

The existing literature about measurements on Bose- 
condensed systems [1,2] concentrates on simple most 
beautiful examples of Bose-condensed states. In this let- 
ter we construct for the first time a measurement theory 
able to predict probability of different density measure- 
ment outcomes on a realistic Bose-Einstein condensate of 
interacting particles evolving under influence of external 
time-dependent potentials. This problem is important 
because in many experiments, like phase imprinting of 
dark solitons [3] or condensate splitting in atom inter- 
ferometers [4] , manipulation of the condensate generates 
substantial dynamical depletion which can qualitatively 
affect measured density patterns. In addition to being re- 
alistic our theory also reveals a beautiful diagonal struc- 
ture hidden in the quantum state of a condensate excited 
from its ground state by a time-dependent potential. 

A generic Bose-Einstein condensate consisting of N 
weakly interacting atoms can be described by a number- 
conserving version of the Bogoliubov theory [5] which 
assumes that most particles occupy the condensate wave 
function. Time-dependent problems can be treated with 
the time-dependent version of the theory where the con- 
densate wave function 0o(i) evolves according to the 
time-dependent Gross-Pitaevskii equation (GPE) and 
the Bogoliubov modes u„i{t) and Vm{t) solve the time- 
dependent Bogoliubov-de Gennes equations (BdGE) 
[5]. In this framework, starting with a condensate in 
the ground state, time-dependent perturbation excites 
the condensate to a state which is formally a time- 
dependent Bogoliubov vacuum [0;,) annihilated by the 
time-dependent quasiparticle annihilation operators 
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Here the operator ag — {(j)o{t)\ijj) annihilates an atom 
in a condensate wave function 0o(t)- The modes Um{t) 
and v^{t) are solutions of the time-dependent BdGE pro- 
jected on the subspace orthogonal to (po, i.e. (Mm|0o) = 
(ti*„|0o) = 0, and normalized so that {um\un) — {vm\vn) — 
Smn- In the time-dependent vacuum |Ofc) the reduced sin- 
gle particle density matrix 
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p{f,r') = Nroif^Mr') + E (2) 

m=l 

has common eigenstates with the operator 

oo 

dp = ^ \Vm){Vm\ = ^diV«|0;)(C| . (3) 
m=l a 

A mode |</)* ) is a non-condensate eigenmode of the den- 
sity matrix (2) occupied on average by dNa particles. 
Thanks to the completeness relation 
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valid in the subspace orthogonal to \4>o), the operator 
dp* = |^'TO)('*^ml common eigenstates with the 

oo 

J2 \Um){Um\ = 1 + dp* = ^(l + dA/-„)|<^a)((Aa|. (5) 

m— 1 a 

Furthermore |^a)'s allow us to get a diagonal form of 

A = Y\Vm){Um\ = J2^-\'>*-){M, (6) 

m a 

which is a map between the subspace orthogonal to \(j)o) 
and the subspace orthogonal to \(Pq). Indeed, the prop- 
erty {um\un) — (vmlvn) = ^nm of the Bogoliubov modes 
and Eq.(4) imply a quasi-commutator dpA = Ad/3* from 
which follows the diagonal form in Eq.(6). What is 
more, owing to dp*{l + dp*) = A*A, one gets |Aa|^ = 
dA^c«(l + dNa). In Ref. [6] we have shown that in the 
stationary case the Bogoliubov vacuum has the following 
form 
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\0b) ~ I al&l + Y Zo.,iajl I |0) , (7) 

a, 0=1 

where a^'s annihilate atoms in single particle states 
orthogonal to the condensate wave function 00 . Matrix 
Za,0 fulfills the equation 
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(8) 



which is equivalent to the set of the conditions bm\Ob) = 
0. In the stationary and time reversal invariant case all 
modes can be chosen as real functions so that the ma- 
trix Zff^a is real symmetric and can be diagonalized by 
an orthogonal transformation [6,7]. In a generic time- 
dependent case, or when a stationary 0o breaks the T- 
invariance like in e.g. vortex state, the matrix Z is com- 
plex and symmetric. However, using only the general 
properties of dp and A it is easy to show that the choice 
of |0a)'s for the modes in Eq.(7) makes the time- 
dependent Bogoliubov vacuum diagonal 
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where the eigenvalues of the matrix Z are 
Afv — — : . 



(9) 



(10) 



Indeed, an action of '^mi^al'^m) on both sides of Eq. (8) 
yields ((/>* |dp|(;ip = ((^* |A|(?i„')^c<',/3 , equivalent to 
Eq.(lO). 

Hence the procedure of time evolution of the Bogoli- 
ubov vacuum state in the particle representation becomes 
extremely simple: i) First one has to evolve Bogoliubov 



modes (m„ 



together with the condensate wave func- 



tion 00 which is easy because different modes evolve inde- 
pendently from each other, ii) Next one has to diagonal- 
ize the operator dp = |'!^m)(^^m| (diagonalization of 
1+dp* — "^jn \'U'm){um\ turns out to be better convergent 
numerically) in order to get eigenvalues dA^^ and eigen- 
vectors \<pa)- iii) Finally one calculates A^ = (0* |A|0„) 
and from Eq. (10) obtains values of Aq,. Note that there 
is a one to one correspondence between moduli of Aq. and 
the eigenvalues of the single particle density matrix dA^„: 
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However, the phases of A's cannot be determined from 
the single particle matrix only, but they also depend on 
the eigenvalues A^ = (0* |A|0a), see Eq.(lO). Finally, 
once the phases are known it is convenient to make the 
transformation 
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In the following we skip the primes. 

For large A'^ the diagonal state (9) can be rewritten as 
a superposition over iV-particle condensates 
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\N : mi)) ■ (13) 



with the normalized condensate wave functions 



0(f|g) = 



(14) 



The gaussian amplitude gives an accurate vacuum state 
when Ea diVa < ^/N . The gaussian square of the am- 
plitude 



p(5)~n^" 



(15) 



is a candidate for the probability distribution for different 
qa- If this choice is correct, then simulation of the density 



2 



measurement amoLmts to choosing a random set of qa 's 
which determines the measured density 

p{r\q) = N\<t>{f\q)\K (16) 

To show the correctness of Eq.(15), we average p{r\q) 
over P{q) and for diV^ <C iV we get the density 

p{f) = 7V|<^o(r-)|2 + 2(1 -A 

averaged over different reahzations of q. This average 
density should be equal to the expectation value of the 
density operator in the state (13) 

{Q,\'^\r)i^{f)\Q,) = 7V|<^o(rl|^ + 5^dArj0„(r-:)|2 (18) 

with dA^Q, = A^/(l — A^). These two densities arc the 
same for the highly occupied modes with diVc ^ 1 or 
Aa Ri 1~. There are discrepancies for poorly occupied 
modes but those modes give negligible contribution to the 
total density. The probability (15) accurately reproduces 
the single particle density matrix. This accuracy can 
be best explained when we look at the overlap between 
condensates 

{N : (l>{r\q)\N : (t>{r\q')) = e" E J?-"?^)' . (19) 

Condensates become orthogonal on the length scale of 1 . 
For the highly occupied modes the fluctuations of qa ~ 
s/ANa are much greater than the width of the overlap 
and for these modes the overlap can be replaced by a 
delta function 5{qa — q'^)- With the delta overlap the 
condensates are orthogonal and the probability P{q) is 
simply given by the Born rule in Eq.(15). 

To summarize, a density measurement, which is a mea- 
surement of all atomic positions, can be approximately 
simulated in two steps. In the first step a condensate 
wave function (14) is chosen randomly from the gaussian 
distribution P{q) ~ exp(— g^/2dA''Q). In the second 
step atomic positions are measured in the chosen con- 
densate. The first step already gives a smoothed density 
distribution A^|(/)(r|g)p on top of which the second step 
only adds statistical fluctuations. In most applications 
one is not interested in the statistical fluctuations but in 
the smoothed density obtained after filtering out the sta- 
tistical noise, compare the histogram and the solid line 
in Fig. 2. Thus for most applications the density mea- 
surement can be very efficiently and accurately simulated 
with only the first step of the procedure which immedi- 
ately gives a smoothed density p{r\q) = N\(f){r\q)\^ with 
randomly chosen g's. In the following we give an exam- 
ple of this procedure: density measurement after phase 
imprinting of a dark soliton. 

To describe evolution of the condensate wave function 
we solve the time-dependent GPE that in harmonic os- 
cillator trap units reads 



1 

idt(l>o = -^dl(t>o + + g\M''^o + Vit, x)cl)o • (20) 

Here V{t, x) is an external potential created by the laser 
beam. As the effective ID interaction strength we choose 
g = 7500 corresponding to parameters of the Hannover 
experiment [3]. The dark soliton is a kink <j)o ~ tanh(a;/^) 
with the wave function winding a phase of tt as a; goes 
from negative to positive. At the same time the den- 
sity drops to zero at x = explaining the name dark 
soliton. Experimentally the soliton is excited by a short 
laser pulse which imprints on the wave function a phase 
that changes by tt when one goes from negative to posi- 
tive X. To simulate the phase imprinting we started with 
the ground state of the stationary GPE for a condensate 
in the harmonic trap, 4>o{x,0), and applied (similarly as 
in Ref. [8] where they use rather small g = 200) the per- 
turbation V(t,x) = Vo[l -I- tanh(x/Zo)] (where Vq = 140, 
lo = 0.05) that lasted for 0.034. Then (j)o{x,t) was 
evolved up to t = 0.15 (or 1.7ms). The density of the 
resulting wave function is plotted in Fig 1(a) where the 
dark soliton is clearly visible in the center of the trap. 

The Bogoliubov modes solve the time-dependent 
BdGE 

idtUm = {Hop + g\4>o{x,t)\'^) Um + 94'o'^m, 

idtVm = - {HgP + glM^^t)]^) - g(l>fum, (21) 

where Hqp = + x^ /2 + V{x,t) + g\(f>o(x,t)\^ . We 

have evolved the Bogoliubov modes, starting with the 
eigenmodes of the stationary BdGE for a condensate in 
the ground state, up to < = 0.15. The single particle 
density p{x,t = 0.15) is plotted in Fig. 1(b). In Fig. 1(c) 
the densities of condensate, noncondensate atoms and the 
total single particle density are plotted. At both values of 
g (the realistic g = 7500 here and the g = 200 in Ref. [8]) 
the single particle density p{x,t) shows that minimum 
at the soliton location is gradually filled with particles 
depleted from the condensate. These dynamical studies 
confirm earlier predictions that a static soliton is going 
to fill up after a few milliseconds [9]. 

After the evolution of the Bogoliubov modes was com- 
pleted we calculated the eigenvalues Aa and the non- 
condensate eigenmodes (j)a{x) at the final time. It turns 
out that the soliton notch is filled up with atoms depleted 
to only one mode: the (f>i (x) with the largest dA^i = 273. 
In Fig. 2(d) the total depletion density Y^„^ \vm{x)\'^ (solid 
line) is plotted together with a barely visible dashed plot 
of the dNi\(j)i{x)\'^ . As the depletion in the soliton notch 
comes only from a = 1 and because we would like to know 
what will be seen in the soliton notch in a single exper- 
iment we truncate the final time-dependent Bogoliubov 

vacuum to ( ajaj -|- Aia|a{ j |0). We made exact sim- 
ulation of density measurement in the truncated vacuum 
using the algorithm of Ref. [1] and Ref. [9]. A typical 
outcome is shown as the histogram in Fig. 2. It turns 
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out that every observed histogram can be very well fit- 
ted with the density p{x\qi) = A^|^o(a;) + qi4ii{x)/\/N\'^ 
where the only free parameter is qi. The histogram in 
Fig. 2 is fitted by the density p{x\14.5) shown in the figure 
as the solid line. This is in agreement with our measure- 
ment theory predicting the smoothed histograms to be 
p{x\qi) with a random qi. 

In conclusion, we derived a convenient diagonal form of 
the time-dependent Bogoliubov vacuum which greatly fa- 
cilitates simulations of quantum measurements on Bose- 
condensed systems. 
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FIG. 1. Panel (a) shows condensate density, i.e. 

N\4>o{x,t)\'^ , at the time t = 0.15. In panel (b) single parti- 
cle density p(x,t) = N\(ptt(x,t)\^ + \vm{x,t)\^ is plotted. 
Panel (c) shows single particle density p(x,t) (solid line), con- 
densate density N\(l){){x, t)\^ (dashed line) and density of non- 
condensate atoms \ vm{x,t)\^ (dotted line) in the vicinity 
of the trap center. In panel (d) ^ \vm{x,t)\^ (solid lino) is 
plotted together with dAfi|0i(x,t)|^ (dashed line barely visi- 
ble behind the solid line) where 4'i{x,t) is the most occupied 
mode (dA^i w 273) orthogonal to (f)o{x,t). The parameters 
correspond to the Hannover experiment where N = 1.5 • 10^ 
(time t = 0.15 corresponds to 1.7 ms). 
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FIG. 2. The histogram is the normalized density distri- 
bution measured in a single exact simulation of the density 
measurement. The solid line is a fit to the histogram with 
one of the predicted p{x\qi) = N\<j>o{x) + qi4>i{x) / \/N\^ with 
gi = 14.5. The dashed line is the total average density p{x) 
which is equal to an average over densities measured in many 
experiments. 
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